Proposal for a Topological Plasmon Spin Rectifier 
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We propose a device in which the spin-polarized AC plasmon mode in the surface state of a 
topological insulator nanostructure induces a static spin accumulation in a resonant, normal metal 
structure coupled to it. Using a finite-difference time-domain model, we simulate this spin-pump 
mechanism with drift, diffusion, relaxation, and precession in a magnetic field. This optically-driven 
system can serve as a DC "spin battery" for spintronic devices. 



The "topological insulator" (TI) is a class of 
strongly spin-orbit-coupled materials with topologically- 
protected, chiral surface states (of opposite group veloc- 
ity and spin) crossing a bulk electronic gap.jlj In these 
surface states, spin and momentum are perfectly related, 
so that the momentum asymmetry caused by an induced 
charge current necessarily spin-polarizes the system. 

Recently, Raghu et ai0 have theoretically investigated 
the properties of resonant electromagnetic excitations of 
charge carriers (plasmons) in the surface states of these 
materials. Because of the large oscillating current in- 
duced at resonance by coupling to incident periodic elec- 
tric fields, a potentially large time-dependent spin density 
is created. It was suggested that these spin-polarized car- 
riers could be used to generate spin currents through non- 
magnetic materials in unconventional "spintronic" de- 
vices, but it is not clear how this could be done since the 
time-average of the instantaneous polarization is zero. 

In this Letter, we propose a means to rectify the 
plasmon-induced AC spin accumulation in TIs by spa- 
tially segregating them in an adjacent resonant plasmonic 
nanostructure. By simulating time-dependent spin drift, 
diffusion, and relaxation, we show that non-equilibrium 
spin polarizations created by the plasmon in the TI can 
be preserved in steady-state at the edges of the nonmag- 
netic metal structure, creating a "spin battery" to poten- 
tially drive spin-electronic devices. 

A schematic illustration of the coupled TI-NM plasmon 
spin rectifier is shown in Fig. 1. Because the perpendicu- 
lar component of current must be zero at any boundary, 
the maximum current (and hence maximum spin density) 
is at the center of the TI structure during plasmon excita- 
tion. At this location, a second plasmonic nanostructure 
with identical resonance, made from a "trivial" nonmag- 
netic (NM) conductor makes contact. Spin generated in 
the TI can thus flow via diffusion across the interface into 
the NM. When spin-up (S y > 0) is generated by a pos- 
itive electron current in the TI, it diffuses into the NM 
and is carried by drift in the +x direction. Half an optical 
cycle later, spin down is generated by a negative electron 
current in the TI and carried by drift in the -x direction 
in the NM. In steady-state, this spatial segregation will 
be maintained in the presence of constant spin diffusion 
and relaxation. 
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FIG. 1: Spin rectification scheme, (a) During one half plas- 
mon cycle, electrons in the chiral topological surface state 
move in the -fx direction and are necessarily in the S y -up spin 
state. They first diffuse across the interface with a plasmon- 
resonant normal metal (1) and then are carried by drift in 
the +X direction (2). (b) During the other half of the plas- 
mon cycle, opposite spin is carried in the opposite direction, 
thus spatially rectifying the time-varying spin density in the 
topological plasmon. 



To substantiate this qualitative description, we numer- 
ically simulate the spin-rectification mechanism by iter- 
ating the one-dimensional spin drift-diffusion-relaxation 
equation in the NM 



dS y 
~~dt 



d 2c ! 



dx 2 dx 



[v(x,t)S y ] - S v /t s 



(1) 



with a Crank-Nicolson technique until steady-state is 
reached. Here, v(x,t) = Vd s'm(irx/L) cos(wi), where Vd 
is the maximum drift velocity of carriers in the NM, L 
is the length of the nanostructure, and u is the plasmon 
radial frequency. D is the spin diffusion coefficient and 
t s is the spin relaxation time. As shown in Fig. 2, zero 
spin-current is maintained at the boundaries of the one- 
dimensional lattice both by the explicit form of this drift 
velocity function (v(0,t) — v(L,t) = 0) but also by mir- 
roring the spin density across the boundary to eliminate 
spin diffusion. Spin injection into the center of the NM 
at x = L/2 due to diffusion from the TI is modeled by a 
source proportional to S^ I cos(ujt) — S y (L/2), where Sy 1 
is the maximum spin density in the TI surface state. All 
parameters for the NM are taken from known values for 
Al at room-temperature: D=40cm 2 /s, T s =65ps.[l] 

We now estimate the plasmon characteristics, treating 
the surface state of the TI dot as an oblate spheroidal 
plasmonic element of diameter L and effective thickness 
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FIG. 2: One-dimensional simulation of spin rectification in 
the NM includes finite-differences Crank-Nicolson solution of 
the spin drift-diffusion-relaxation equation [Eq. {JJ] subject 
to zero spin-current boundary conditions and diffusive injec- 
tion of spin from the TI plasmon into the center of the NM 
plasmon. 



t along z. The electric polarization of the plasmonic ele- 
ment can be expressed asQ 
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where Eq is the incident field, n x is the depolarization 
factor of the spheroid for Eq, xq is the spatial charge 
carrier displacement and ri2 is the 2-dimensional density 
of the surface state. The effective static dielectric exter- 
nal to the surface state e we take as the average between 
vacuum and the bulk TI ss50. The effective dielectric 
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where 



of the surface state e is taken as e , , 

-^f. For t « L, n»« ft.® 
The electric polarization is maximum when the denom- 
inator of Eq. @ is minimized, which gives the plasma 
resonance frequency uj 2 — ™ L ^i . [2| At 
then given by 
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FIG. 3: Spatial spin density evolution in a 1-d model of 
a normal-metal plasmonic structure coupled to a TI in res- 
onance at 700GHz. (a) shows the cycle-average spin den- 
sity and at several times throughout the cycle with period 
T p ~ 1.4ps; (b) shows the full evolution with time on the 
vertical axis. Throughout a full cycle of plasmon oscillation, 
a steady-state spin accumulation of equal and opposite sign 
is apparent, thereby constituting DC spin rectification of the 
AC spin generation from the topological plasmon. 
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where is the momentum scattering time; we assume 
LJTk < 1. The maximum drift velocity is simply Vd — ujxq. 

For L=l /im and Ep=100 meV, the resonant frequency 
/ = ui/2ir «700 GHz. A relative amplitude of xq/L — 
0.1 requires an incident peak power of 100 W focused 
to 1 mm, which is well within the range of pulsed THz 
radiation sources. Under these conditions, Vd is on the 
order of 10 7 cm/s; we assume the modes in the TI and NM 
are electromagnetically equivalent and use this value in 
the following simulations. 

In Fig. 3(a) and (b) we show the simulation results of 
S y (x) during one complete plasmon cycle using the pa- 
rameters given above after approximately 2.5ns simula- 
tion time. These simulations use a spatial discretization 
Ax = lOnm and time step of At = 10 fs, such that the di- 
mensionless parameters DAt/Ax 2 , VdAt/Ax < 1, which 
is needed for numerical stability. Although the spin dif- 
fusion into the center of the NM plasmonic structure is 
harmonically oscillating in time, the spatial symmetry 
breaking provided by the plasmon-induced drift gives rise 



to a static spin accumulation which diffuses to the edges 
of the NM. The cycle-averaged spin density clearly shows 
opposite values of S y on either side, constituting a polar- 
ization splitting and spin-battery "potential" AS y — 0.231 
relative to the maximum spin density in the TI, Sy 1 . 

It is useful to determine the effect of changing param- 
eters on this polarization splitting. In Figs. 4(a)-(e), we 
plot ASy in black as a function of Vd, T, L, D n , and 
T p = 2ir/uj around the fixed values given above. Clearly, 
to maximize the spin battery output one desires intense 
electromagnetic illumination causing large drift veloci- 
ties, a large spin lifetime in the NM, and low plasmon 
frequencies. 

The relaxation-time approximation can give a heuris- 
tic prediction of all the trends shown in Figs. 4(a)-(e). 
Spin-up and spin-down regions are separated by periodic 
carrier drift by a distance of w VdT p . The loss due to 
(?) spin relaxation in the region towards the boundary 
at a rate s» -^f-{Lj2 — VdT p /A) and (ii) spin diffusion 

AS 1 

from spin- up to spin-down regions equal to D — =f must 
be balanced in steady state by spin injection. The lat- 
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consequences of lower aspect ratios can be simu- 
lated with a two-dimensional extension of the spin 
drift-diffusion model described above. 

2. The spatial overlap between TI and NM nanostruc- 
tures in our simulations is modeled via spin diffu- 
sion into a single finite-difference element. In prac- 
tice, the spatial overlap should be smaller than the 
carrier drift distance Vd/f (here 143nm, well within 
the range of electron-beam lithography). 

3. Matching the plasmon resonance frequencies be- 
tween the TI and NM structures is crucial, and 
is aided by (anti-)symmetric mode coupling and 
phase-locking. (For small structures where L « 
c/lo, only the symmetric mode is excited by the in- 
cident field.) In practice, the TI resonance may be 
tuned by electrostatic gating if the gate itself does 
not prohibitively interfere with the plasmon mode. 
This gate can also be used to adjust Ep below the 
Dirac point where there will be a second resonance 
(of spin-polarized holes) . 




FIG. 4: Dependence of spin accumulation AS y on trans- 
port parameters in the NM plasmon nanostructure: (a) drift 
velocity Vd, (b) spin lifetime t s , (c) length L, (d) Diffusion 
coefficient D, and (e) plasmon period T p . Black lines are the 
results of Crank-Nicolson simulations of Eq. ((T]) and grey 
lines are the heuristic predictions of Eq. (|4)| . Vertical dashed 
lines indicate values used in calculations shown in Fig. 3. In 
(f), a perpendicular magnetic field B z causes spin precession 
and dephasing. 

ter proceeds by diffusion and hence is proportional to 
(Sy 1 — ASy/2). Therefore, in the weak injection limit, 
we expect that the spin-battery strength scales as 

iS »=< ^ r + W, [-2-^)1 ■ <4) 

Favorable comparison of the trends predicted with this 
method is made to the numerical model in Figs. 4(a)-(e) 
(grey lines). 

Our simple spin drift-diffusion-relaxation model sug- 
gests some caveats. For instance: 

I. The validity of a one dimensional approximation re- 
quires the diffusion time across the interface should 
be faster than half the plasmon cycle. For the pa- 
rameters used here (D=40 cm 2 /s and f=700GHz), 
the allowed length yj D / f= 38nm. This will require 
nanostructures of high aspect ratio. The detailed 



4. The effect described here is dependent on the 
strength of the oscillating spin density Sy 1 in the 
TI driving diffusion into the NM. The spin polar- 
ization P w AE/Ep = eEyX/Ep = Vd/vF (where 
AE is the Fermi-level asymmetry between spin-up 
and -down states, A is the mean-free-path, and vf 
is the Fermi velocity) is « O.f here. Significant 
increases in this value may be limited by energy 
relaxation through strong optical phonon coupling; 
for instance, in Bi2Se3 at AE w8 meV.jBJ 

To experimentally confirm the presence of the ex- 
pected spin accumulation in the NM, we propose to 
detect nonzero AS y with ferromagnetic tunnel voltage 
probes. @ To rule out spurious signals due to asymmet- 
ric tunnel barriers and current rectification/photovoltaic 
effects across the contacts, a magnetic field B z perpen- 
dicular to the surface (in the z direction) will be used to 
induce spin precession and subsequent dephasing (Hanle 
effect), suppressing this voltage signal. [6( By modeling 
evolution of both S y and S x spin components, and by in- 
corporating spin precession in Eq. (1) by adding a term 
_ m^L> g x (S x x + S v y),$ (where g = 2 is the electron 
spin g- factor, fj,s is the Bohr magneton, and Ti is the re- 
duced Planck constant) we have simulated this suppres- 
sion, shown in Fig. 4(f) for both AS y and the spin signal 
in a perpendicular, in-plane direction AS X . In general, 
longer spin lifetimes will result in stronger low-B-field 
suppression. 

Finally, we wish to point out that this plasmon-induced 
spin generation effect is not necessarily confined to TIs. 
For example, the spin-Hall effect has been observed in 
extrinsically doped n-GaAs.0] Excitation of plasmons in 
GaAs will cause a transverse AC spin accumulation; sub- 
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sequent diffusion of these spins into resonant nanostruc- 
tures with limited spin-orbit interaction strength is ex- 
pected to result in spatial rectification as in the above 
scenario. 
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